HD-THEP-99-26 



The field strength correlator from QCD sum rules 
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The gauge invariant two-point correlator for the gluon field strength tensor is analysed by means of the QCD 
sum rule method. To this end, we make use of a relation of this correlator to a two-point function for a quark- 
gluon hybrid in the limit of the quark mass going to infinity. From the sum rules a relation between the gluon 
correlation length and the gluon condensate is obtained. 



1. Introduction 

An important role in non-perturbative ap- 
proaches to QCD is played by the gauge invariant 
non-local gluon field strength correlator 



b 
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(1) 



where the field strength F^^ = d^A^ - d^A^;^ + 
gf^bcj^b^j^c^ z = 2/ — a; and V denotes path order- 
ing of the exponential. In general, the gauge in- 
variant field strength correlator could be defined 
with an arbitrary gauge string connecting the end 
points X and y, but in this work we shall restrict 
ourselves to a straight line. 

It is the basic ingredient in the stochastic model 
of the QCD vacuum and in the description 
of high energy hadron-hadron scattering H-|6) . In 
the spectrum of heavy quark bound states it gov- 
erns the effect of the gluon condensate on the level 
splittings ||7|-|l0|] and it is useful for the determi- 
nation of the spin dependent parts in the heavy 
quark potential • 

The correlator can be related to a correlator of 
a colour singlet current composed of a (fictitious) 
infinitely heavy octet quark and the gluon field 



*Invited talk at the Euroconference on Quantum Chromo- 
dynamics (QCD'99), Montpellier, July 1999 



strength tensor. This fact has already been em- 
ployed in ref. |Q in order to calculate the per- 
turbative corrections by means of Heavy Quark 
Effective Theory (HQET) . In this work we again 
use this relation to estimate the correlation length 
from QCD sum rules using as ingredients the 
value of the gluon condensate and the results for 
the perturbative calculation. 

2. The field strength correlator 

Instead of dealing with the non-local string op- 
erator which makes a calculation rather tedious 
one can replace the string by local heavy quark 
fields and use the methods developed in HQET. 
To this aim, we introduce an infinitely heavy 
quark field in the octett representation, h°'{x), 
which is constructed from the field analo- 
gous to HQET by 



h''{x)= lim i(l+^)e 



imnvx r^a 



(2) 



V is the four-velocity of the heavy quark. The 
propagator of the free heavy quark field in coor- 
dinate space is given by 



S{z) 



{Q\Tmy)h\x)m 



(3) 



where v'^ is the zero-component of the velocity. 
The 5-function constrains the heavy quark on a 
straight hne. With the effective HQET action 
Seff = J dx hiv^Di^^h , the following equation 
can be shown analytically JTst : 

{0\T{h''{y)h\x) e'^-'^^ llO) 

= S{z) (0|7'e^-^°""^X^*^^(^+*^)|0) . (4) 

The physical picture of this result is a heavy 
quark moving from point x to y with a four- 
velocity V, acquiring a phase proportional to the 
path-ordered exponential. The limit of mg — > oo 
is necessary in order to constrain the heavy quark 
on a straight line and in order to decouple the spin 
interactions which are suppressed by a power of 
1/niQ and can therefore be neglected. By intro- 
ducing a new correlator 'D(z), we can now express 
our correlator (^) in terms of heavy quark fields 

^fii^pa (^) 



S{z) Dfj,i,pa{z) 



(5) 



which establishes the relation between the field 
strength correlator and HQET. 

3. The sum rules 

Our next aim is to evaluate this correlator in 
the framework of QCD sum rules jlj] and in that 
way obtain information on the correlation length 
of the field strength correlator. We may view the 
composite operator [g gh"" F1^^){x) as an interpo- 
lating field of colourless quark gluon hybrids and 
evaluate 'D^^p„{z) by introducing these as inter- 
mediate states in the absorption part of 'Dp.upaiz)- 
The lowest lying state will govern the long-range 
behaviour and hence the inverse of its energy is 
the correlation length. 

For the sum rule analysis it is preferable to 
work with the correlator in momentum space. 



Thus we define 



2?p.p.H = ^ / dze^'^%Q\T{glF-Mh''{y) 



xFl[x)h\x)}\Q), 



(6) 



where w = vq is the residual heavy quark mo- 
mentum. 

For a sum rule analysis the states have to be 
classified according to different quantum num- 
bers. The projections onto the two independent 
subspaces is done with 



V~{w) EE g^'Pv''v''V 



^upa ('^) 

(5^V"-2g^^t;^i,-)P^,p.H, (7) 



where 'D~ contains a vector and T>^ an axialvec- 
tor intermediate state. 

We model the correlators by a contribution 
from the lowest lying resonance plus the pertur- 
bative continuum above a threshold sq- Inserting 
the matrix elements and performing the heavy 
quark phase space integrals one obtains 



,-T I fT|2 



X — w — ie ' 



(8) 



where E represents the energy of the glue around 
the heavy quark, are the hadronic matrix el- 
ements and are constants. The spectral den- 
sities are defined by p^(A) = l/7rImI>^(A -|- ie) 
and are known at the next-to-leading order . 

After Fourier transformation to coordinate 
space the above representation reads: 



(9) 



Since the heavy quark propagator factorises, we 
identify the expression inside the brackets with 
our original correlator 'D{z). The long-range be- 
haviour will be dominated by the term containing 



the exponential with the energy E. Therefore the 
correlator decays exponentially and the correla- 
tion length is given by 1/E. 

Now we turn to the theoretical side of the sum 
rules which is obtained by calculating the corre- 
lator of eq. (H) in the framework of the operator 
product expansion ||l^ , |l6|] . 

The perturbative contributions in momentum 
space have the form 



+aip^o+P2iL+P^2L')\, (10) 

where a = as/Tr, L = ln(— 2w//i) and the coeffi- 
cients pfj can be found in rcf. [0 . 

Essential for the sum rule analysis arc the con- 
tributions coming from the condensates. In our 
case the dimension three condensate (hh) van- 
ishes since the quark mass is infinite. The lowest 
nonvanishing term is the gluon condensate of di- 
mension four: 

The next condensate contribution would be of di- 
mension six, but we shall neglect all higher con- 
densate contributions in this work and restrict 
ourselves to the gluon condensate. 

The correlators satisfy homogeneous renormal- 
isation group equations. Thus we can improve the 
perturbative expressions by resumming the loga- 
rithmic contributions. Calculation of the first co- 
efficients 7^ for the anomalous dimensions from 
eq. (pel) leads to 



V~p{w) - \ V+p{w) = - ^ {aFF) . (11) 
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0. 
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(12) 



The correlator P~(w) which corresponds to the 
vector intermediate state does not depend on the 
renormalisation scale at this order. 

In order to suppress contributions in the disper- 
sion integral coming from higher exited states and 



from higher dimensional condensates, it is conve- 
nient to apply a Borel transformation Bt with T 
being the Borel variable. After renormalisation 
group improvement and Borel transformation all 
the ingredients needed for a sum rule analysis are 
known. Explicit formulas for the expressions can 
be found in ref. (it) . Now we turn to the numer- 
ical analysis. 

4. Numerical analysis 

Let us denote by x(r, Sq) the Borel transformed 
expression for the continuum part xC^: ^o) = 
Bt d\ p{X) / {X — Lu — ie) . After equating the 
phenomenological and the theoretical part we end 
up with the sum rule 

= Vfp+V^^iT)~x^(T,so). (13) 

By taking the logarithmic derivative we get an 
equation for the energy E. 

The analysis shows that the different sign of 
the perturbative and non-perturbative term in 
the 1~ state leads to a stabilisation for the energy 
sum rule, whereas the equal sign in the 1+ state 
destabilises. 

As our input parameters for the case of three 
light quark flavours we use < aFF >= 0.024 ± 
0.012 GeV^ A3// = 0.323 GeV and = 2 GeV. 
To estimate the errors we have varied the scale 
/i, the continuum threshold sq and the gluon con- 
densate. In Fig. 1 we have displayed the energy 
E~ as a function of the Borel parameter T for 
different values of fj, and Sq. A good balance be- 
tween stability for the energy and sensitivity for 
the resonance is found around E^ = 1.5 GeV. 

Including the errors we get for the energy and 
correlation length 



E. 



-= 1.5 ± 0.4 GeV, fl-^; = O.lSI^ 'oa fm. (14) 



In a world without light quarks, i.e. nf — 0, 
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Figure 1. The energy as a function of the Borel- 
parameter T. Dashed curves ^ = 1 GeV: low- 
est So = 1.3 GeV, middle so = 1.5 GeV, upper 
So = 1.7 GeV. Sohd curves ^ = 2 GeV: lowest so = 
1.5 GeV, middle sq = 1.7 GeV, upper so = 1.9 GeV. 
Dotted curves /i — 4 GeV: lowest sq = 1.7 GeV, 
middle so = 1.9GeV, upper so = 2.1GeV. 

the main influence is the expected change of the 
gluon condensate which might increase by a factor 
two to three If we perform an analysis as 

above, we get for Aq/j = 0.250 GeV, < aFF >= 
0.048 ± 0.024 GeV* and so = 2.3 GeV an energy 
and correlation length of 

E-^i = 1.9 ± 0.5 GeV , a^^i = 0.1lt|]:J]^ fm. (15) 

For £'+, the energy of the axial vector, we obtain 
no stable sum rule. The expressions for E~ and 
E^ are equal in lowest order perturbation the- 
ory, higher order perturbative contributions and 
the gluon condensate lead to a splitting in such 
a way that for the same values of so and T the 
expression for E^ is higher than that for 

5. Summary and conclusions 

The analysis of the gauge invariant gluon field 
strength correlator by QCD sum rules allows to 
establish a relation between the gluon condensate 
and the correlation length. In order to apply the 



sum rule technique which consists in the compar- 
ison of a phenomenological ansatz with a theo- 
retical expression obtained from operator prod- 
uct expansion we interpret the gluon correlator as 
the correlator of two colour neutral hybrid states 
composed of a (fictitious) heavy quark transform- 
ing under the adjoint representation and a gluon 
field. The former serves as the source for the 
gauge string in the correlator. 

The value of the lowest intermediate 1^ state 
(the inverse correlation length of the correlator) 
with three flavours can be determined to E^j^ = 
~ 1.5 ± 0.4 GeV and with zero flavours 
to Eqji = l/flQ^; « 1.9 ± 0.5 GeV. The main 
sources of uncertainty are the choice of the con- 
tinuum threshold sq and the value of the gluon 
condensate. 

Though we find no stable sum rule for the axial 
vector state we have from the difference of the ex- 
pressions for the 1~ and 1^ state strong evidence 
for the counterintuitive result that the 1+ state is 
lighter than the vector state. 

The field strength correlator has been calcu- 
lated on the lattice using the cooling technique 
[|l^,|l^ and field insertions into a Wilson loop j2^ . 
For a discussion see ref. 0. Recently m an 
analysis of heavy quarkonium in the framework 
of NRQCD ||2^] a splitting between the vector 
and the axialvector part has been observed in the 
same direction as proposed by the sum rules. 

The sum rule analysis shows that the state in- 
vestigated here namely a gluon confined by an 
octet source has a much higher energy than a cor- 
responding state in HQET. A similar analysis of 
a light quark bound by a source in the fundamen- 
tal representation yielded an energy which is 
by a factor 2 to 4 smaller. This is to be expected 
from general grounds since the case treated 
here is nearer to a glueball than to a heavy meson. 
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